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We consider the visco-elastic response of the electronic degrees of freedom in 2D and 3D topological 
insulators (TI). Our primary focus is on the 2D Chern insulator which exhibits a bulk dissipationless 
viscosity analogous to the quantum Hall viscosity predicted in integer and fractional quantum Hall 
states. We show that the dissipationless viscosity is the response of a TI to torsional deformations 
of the underlying lattice geometry. The visco-elastic response also indicates that crystal dislocations 
in Chern insulators will carry momentum density. We briefly discuss generalizations to 3D which 
imply that time-reversal invariant TI's will exhibit a quantum Hall viscosity on their surfaces. 

PACS numbers: 



A striking feature of a topological insulator (TI) is 
its topological response. The paradigmatic example is 
the time-reversal breaking integer quantum Hall effect 
(IQHE) in 2D which exhibits a Hall conductance that is 
an integer multiple of e 2 //ipQ . More recently it was shown 
that there exist related states in 3D which are time- 
reversal invariant [21 [3] and exhibit a topological magneto- 
electric effect [3] (TME). While the electro-magnetic re- 
sponse of topological insulators is the most well-known, 
in this Letter we consider the visco-elastic response of the 
electronic degrees of freedom in TI's. Namely, we want 
to consider the stress response 

{T^)=K» kt u u + rf iU uu (1) 

where T 4 -? is the stress-tensor, A, r\ are the elasticity and 
viscosity tensors respectively and Uij is the strain tensor. 
Here we show there is a dissipationless viscosity response 
in the topological Chern insulatorjS] (CI) state analogous 
to that found in the IQHE and fractional QHE states jBJ- 
[TU] . While viscosities are normally associated with fric- 
tional dissipation, this viscosity, present only when time- 
reversal symmetry is broken, implies a force perpendicu- 
lar to the fluid motion similar to the Lorentz force. 

In a condensed matter system the electronic stress 
response can be calculated by coupling the electronic 
Hamiltonian to perturbations of the background lattice 
geometry. The topological responses due to geometric 
curvature have been studied in Refs. I12j in the lan- 
guage of quantum field theory anomalies. Alternatively, 
we consider the response of topological insulators to an 
external torsion field. A heuristic understanding of the 
difference between curvature and torsion is that when an 
object traverses a small loop in real space it is rotated if 
there is non-zero curvature, and translated if there is non- 
zero torsion. A familiar manifestation of torsion is a crys- 
tal dislocation. These line-defects are singular sources of 
torsion, analogous to a localized magnetic flux line. For 
example, while dragging an electron around a magnetic 
flux line its wavefunction is multiplied by a ?7(1) phase, 
while for a dislocation line it is multiplied by a transla- 
tion operator along the Burgers vector. The visco-elastic 



response is derived as a linear response to torsional per- 
turbations of the underlying material geometry. We find 
that this response also leads to momentum density local- 
ized on crystal dislocations and we also briefly mention 
the visco-elastic response of 3D time-reversal invariant 
topological insulators (3DTI). 

To understand the torsion response we will often draw 
comparisons to the well-known electromagnetic responses 
of topological insulators which we briefly review now. 
The responses of the CI and 3DTI to external electro- 
magnetic fields are encapsulated in topological effective 
actions i.e. free-energy functionals derived from calculat- 
ing a partition function in the presence of external fields. 
The QHE and TME are encoded in the effective actions 

n p 2 /* 

SefflA,} = — J d 3 xe^A,d„A p (2) 
SefflA,} = ^/ otxBe^F^ (3) 

respectively which are derived from the responses of the 
topological insulators to an external field A p (note that 
n is an integer). The nominal current response is (j 11 ) = 
SSeff/SA^ which gives the QHE and TME when acting 
on Eq. ([2]) and ^ respectively. All known topological 
electro-magnetic responses in various dimensions can be 
described by similar topological effective actions [I]. 

Our primary interest is the 2D CI for which we will use 
a continuum massive Dirac Hamiltonian as a model. To 
calculate the visco-elastic response we couple the mas- 
sive Dirac Hamiltonian to geometric perturbations. Be- 
cause of its spinor nature, the Dirac Hamiltonian does 
not couple to geometry through the metric tensor, but 
instead via the orthonormal triad e a and its inverse e a 
(frame field) and the spin connection Lu a i>. The latin in- 
dex a labels the particular vector of the frame which, 
when expanded in terms of a local coordinate basis 
d/dx 11 — (d t ,d x ,d y ), has components e%. In a lattice ver- 
sion of the theory, the frame is defined by the local orbital 
orientation. The stress response can be thought of as a 
functional of e a and oj a b, but we should not take them to 
be related to each other as they would be in Riemannian 
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geometry [T3]. In particular, we will find that the dissipa- 
tionless viscosity response of the Dirac model is related 
to torsion. In the context of condensed matter physics, it 
is convenient, in fact, to set the spin connection to zero 
such that the torsion is contained in the properties of the 
triad alone. 

The action and Hamiltonian for continuum 2D massive 
Dirac fermions coupled to a frame field are 

S = J d 3 xdct(e)$( Ptl e^ a - m) i/> 
H = Px elT a +p y elT a + mT° (4) 



with a 



0,1,2, r 



-ia x ) and T a 



(a z , tj x , a y ). If the frame field is position independent 
the energy spectrum is simply E± = ±^/p\ + p\ + m 2 
with p a = e l a Pi. This is a gapped insulator when m/0. 
Now we will calculate the off-diagonal response of the 
stress-energy current (analogous to a xy ) due to a pertur- 
bation of the triad e^(x) = 5® + Se^(x) around the trivial 
background. We will see later that the triad has a sim- 
ple interpretation in terms of elasticity theory and pro- 
vides a natural geometric deformation. The stress-energy 
current that couples to the triad is T£ = det ^ e ^ j^- = 

ipp^jaip. We wish to integrate out the massive fermions 
to get an effective action which is a functional of the 
triad. We are only interested in the terms which lead to 
dissipationless transport and we find, at leading order, 



(T a ^)(q) = ^- Vab e^q a I T (m) 



J T (m) 



(5) 
(6) 



where rj a b = diag[l, —1, —1] is the flat-space Lorentz met- 
ric, q is the external momentum, and y = p 2 where p is an 
internal loop momentum. If we Fourier transform back 
to real-space this kernel leads to an effective action 
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I T (m) 



d 3 x e^e a u d v e b p r]a b . 



(7) 



which is similar to Eq. Q, i.e., a Chern-Simons (CS) 
term for the triad. Restoring the spin connection, the 
integral in Eq.|7]) is the Lorentz invariant integral J e a A 
T b r] a b, with T a the torsion 2-form. For reasons we will 
see below, we call this a quantum Hall viscosity response. 

When probed by an electric field the 2D continuum 
Dirac model is notorious for having a half-integer QHE 
(<?xy — sgn(m)e 2 /2/i) which is connected to the parity 
anomaly [T3]. However, when properly regularized, (e.g. 
on a lattice) a xy becomes quantized in integer units, as it 
must for a non-interacting system [5]. In the present case, 
in the continuum limit, the coefficient It(tti) is divergent. 
If we simply cut off the momentum integral at a UV scale 
A then we find It(™, A) = — mA + 2m 2 sgn m + 0(1/ A). 
Comparing to the quantized Hall conductance, this is 



quite different, although from symmetry and dimensional 
analysis there is no choice: this term must break time- 
reversal and thus is an odd function of m. Additionally 
since is dimensionless (unlike A^) this coefficient must 
have units of 1/ [Length] 2 . Hence the leading term is pro- 
portional to m and the only other scale A. The other 
unusual thing is that this term is continuous at m = - 0, 
unlike the Hall-conductance, which jumps. To get physi- 
cally sensible answers for the Hall viscosity, and the Hall 
conductance, which cannot be half-integer, we must more 
carefully regulate the theory. Here, we describe the stan- 
dard Pauli-Villars technique with a set of N massive reg- 
ulator fields, which is appropriate since it preserves all 
the symmetries of the Hamiltonian. The i-th regulator 
field has mass Mj and wave-function renormalization Q 
and we take Mo = m, Co = 1. The regulated Hall con- 
ductance and viscosity are 
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respectively. We can rewrite 

M r °° 



(8) 



I T (M) = - lim ^= I dt t 1 ' 2 \ dy ye~ t( -y +M > (9) 
^0 y/TT J e Jo 

which yields I T (M) = -2M/^e + 2\M\ 2 sgnM + 0(</e). 
We have three physical constraints that will give proper 
renormalization conditions (i) a xy = ne 2 /h, n £ Z (ii) 
Q reg is finite and (iii) if a xy = then £ reg = 0. We know 
that when the Dirac mass switches sign we go through a 
phase transition between a trivial insulator and a topo- 
logical insulator. The sign of m that gives the topolog- 
ical insulator is regularization dependent, and without 
loss of generality we pick m > to be the non-trivial 
insulator. This means that for m < we require both 
&xy = Creg = 0, which is the origin of the third condi- 
tion. A solution for Mi,Ci under the constraints above 
is always possible. In all cases, one finds that in the 
non-trivial phase we have 
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If we had chosen m < to be the topological phase then 
the signs of both coefficients would be flipped. Note that 
we have restored the units in the viscosity response. The 
dimensions of this coefficient are angular momentum den- 
sity, which is equivalent to momentum per unit length, 
and the units of dynamic viscosity (force/ velocity). Com- 
paring to the value of the Hall viscosity for the IQHE[6l[7] 
Ciqhe — 8 J^2 we see a similar structure coming from the 
length scale endowed by the time-reversal breaking field. 
Also, the viscosity is continuous in the limit m — > anal- 
ogous to the B — > (l B — ► oo) limit of the IQHE. 

To make contact with the literature on the IQH vis- 
cosity we will re-derive the Hall viscosity for the CI state 



using an adiabatic transport calculation. This can be 
carried out by putting the Dirac equation on a torus and 
calculating the adiabatic curvature due to shear deforma- 
tions (equivalently deformations of the modular parame- 
ter t) of the torus [5] . Consider a square torus, made in M 2 
by identifications (x, y) ~ (x+a, y+b) with a,b € Z, with 
fixed unit volume, and consider area preserving diffco- 
morphisms, corresponding to spatial metrics of the form 
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(11) 



The basis vectors and the spatial part of the triad are 
ei = y/ndx, e 2 = -= (-nd x + d y ) (12) 

V T 2 



(dx - Tidy), e 2 = y^dy (13) 



respectively, and the Hamiltonian is 

m V 



H 



V -m 



(14) 



where V = ^= (f Pl - p 2 ) . We define VV = \ \V\\ 2 . 

We consider the ground state in which all of the neg- 
ative energy states ip-(pi,p2]T) are occupied. The adia- 
batic connection can be calculated from the explicit form 
of the single-particle wavefunctions and we find 

A = i ip_(m,n;T)dip^(m,n;T) 

nGZ 

E f(\\n 2 )U\n(Z) (1 5 ) 



m,n£Z 



where 



f(\\n 2 ) 



(m 2 + ||7'|| 2 ) 1 /2 



(16) 



and where the sums are over the discrete quantized mo- 
menta on the torus. This gives the adiabatic curvature 



p = .drA^^ p2jlmn 



If we convert the sum into an integral we find 



(17) 



F = 



,dr A dr It(w) 



which yields the same (unregulated) viscosity as above. 

We will now give a physical interpretation in terms 
of conventional elasticity fields [15]. Assuming we have 
an elastic medium, we can pick a reference un-displaced 
state and define a (space-time) displacement field u a {x). 
Then the triad can be written as e' 



5 a 



«)° where 
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FIG. 1: (a) Chern insulator deformed by a dislocation - 
antidislocation pair, separated in the y-direction. For each 
dislocation, the momentum-density is in the direction of the 
Burgers vector, (b) Chern insulator on a cylinder with a 
(non-uniform) dislocation threading the cylinder. Local dis- 
placements are shown by red arrows. This gives rise to a 
momentum-current response along the cylinder which carries 
a momentum component parallel to the Burgers vector of the 
threaded dislocation, i.e., parallel to the red arrows. 



we assume that u° = i.e. e 

velocity field w A — v A and the spatial distortion tensor 



dt. Now «j° contains the 



,,A 



w\* where A = 1, 2 and i = x,y. This formulation of the 
triad in terms of the distortion tensor is consistent with 
the usual definition as can be seen by calculating the spa- 



tial metric — efef&AB — Sij + wtj + Wji + wfw 
which matches the metric from elasticity theory [Fo] 
stress-energy response from Eq. ([7| is 



"A 

The 



B 



(19) 



Since e° does not enter, this simplifies to a momentum- 
density (T^) = CregVAB^dief & n d a momentum-current 
(T A ) — C,regTlABf lvp d v e^ . These satisfy the continuity 
equations d t (T^) = —di{T\). Restricting ourselves to lin- 
ear elasticity theory we can freely switch between frame 
(a) and local coordinate (/i) indices in the response equa- 
tions. Also, space-time indices are raised/lowered using 
the unperturbed metric. 

For (T Q °) ^ 0, u A cannot be single-valued: it is a dislo- 
cation with Burgers vector b A at Xo, for which e lJ diiu A = 
e^dief — —b A S^(x~K.a). Thus, the momentum-density 
(18) response simplifies to (T oj ) = -( reg J2 m & j m) ^ (2) ( x ^ x »Ti) 
where the x m are the locations of dislocations and b^ 



is the Burgers vector of the m-th dislocation. For a lat- 
tice system, the dislocation is the fundamental quantized 
unit of torsion since transporting an electron around a 
defect translates the wavefunction by a multiple of a lat- 
tice constant. An edge dislocation with |b| = cfj contains 
a momentum density of =± along the direction of b. 



w" 



du a /dxp is the distortion tensor [T5]. To simplify As illustrated in Fig. [lja), this response is a momentum 
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density bound to a torsion 'flux' analogous to charge den- 
sity bound to an electromagnetic flux in the bulk of a CI. 
Note that Fig [I] is heuristic, since a realistic edge dislo- 
cation is not simply a cut into the material. 

The physical interpretation of the momentum current 
response (Ty) is not as simple because it is more diffi- 
cult to picture a torsion electric field. In the 2D plane, a 
moving dislocation (torsion flux) will generate a torsion 
electric field via the analog Faraday effect. Since we have 
seen that dislocations naturally carry a momentum den- 
sity, moving it will generate a momentum- current density 
as per the response equation. In fact, the momentum- 
current due to the moving dislocation is being carried 
perpendicular to the induced torsion electric field. 

Another realization of the momentum current response 
is obtained by using another instance of the Faraday ef- 
fect: roll the CI into a cylinder and then insert a torsion 
flux into the cylindrical hole. This can be thought of 
as threading a dislocation into the hole of the cylinder 
so that any particles traveling around the hole will be 
translated by the Burgers vector b a of the threaded dis- 
location. Changing b a as a function of time creates a 
torsion electric field the same way that a changing mag- 
netic flux causes a circulating electric field. There is one 
key difference with the electromagnetic case: in order to 
preserve area the threaded torsion flux cannot be uni- 
form and the translation must average to zero over the 
length of the cylinder as illustrated in Fig. [ijb). Thus, 
the natural experiment is a torque experiment where a 
cylinder of CI is twisted. This is equivalent to threading 
a dislocation with a position dependent Burgers vector. 

The formalism developed here is a natural generaliza- 
tion of classical elasticity theory. If w® = on the bound- 
ary, we can rewrite the effective Lagrangian (Eq. (j7|)) as 



C 



, ff = ^l e ^P VabW a duW b 



^e^V A {u^d v u p x + M ua 8 u M pX + 2M tl „d v u pX ) 



1 / du^ du v 

2 \dx u dxn 



dx. 



du v 
dx u 



within linear elasticity theory; and M M „ are the strain 
and rotation tensors respectively. The first term is the 
torsional viscosity which is the Lorentz invariant version 
of the QH viscosity [BJ. The stress-energy tensor response 
(T/iv) is not necessarily symmetric and thus does not 
fit in classical elasticity (which would be independent of 
Mpv). It is natural to interpret the stress-response within 
micropolar (Cosserat) elasticity theory which takes the 
local rotational degrees of freedom of the medium into 
account [TBJ ET] . The distinction here is clear since Dirac 
fermions couple directly to the triad and not the metric, 
and the spinor nature of the Dirac equation gives local 
rotational degrees of freedom to which the triad couples. 

Finally, we briefly mention two interesting 3D general- 
izations, the details of which will be presented elsewhere. 



The first is an anisotropic extension to 3D with the form 



Seff K] _ } 



d 4 xe^e a u d p e b a r, ab (20) 



where £ M is a vector of viscosity coefficients which is anal- 
ogous to the 3D IQHE[TH]. IQHE or QAHE states which 
are 'stacked' along a direction perpendicular to the vec- 
tor exhibit the viscosity response in Eq. (20 1. This 



action is basically equivalent to the 3D viscosity response 
of Ref . (TO] . For a 3D strong TI we find 



S, 



1 



effl 



d 4 x ^e^d^dpel^. (21) 



which is a total derivative unless is not a constant. 
Hence, on the surface of a 3DTI (where Q 3 ' has a domain- 
wall) there will be a dissipationless visco-elastic response. 
This is expected since the surface also contains a QHE. In 
gravity theories with torsion the isotropic term is known 
as the Nieh-Yan term|19j and the anisotropic term is the 
(anisotropic) extension of the triad CS term. We leave 
open the question on how to experimentally measure this 
response in 2DEGs or TIs. Unlike electric charge, mo- 
mentum is not conserved when translation symmetry is 
broken, as in any realistic material. Thus questions about 
possible quantization or the nature of the topological re- 
sponse are still unanswered although interesting discus- 
sions of physical implications have been given [9l [TO], 
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